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1. Introduction 

This talk critically reviews recent lattice QCD results relevant for kaon phenomenology, as 
well as the methodology that is used to obtain them. The focus is on full QCD calculations, which 
account for the effects of light sea quarks either partially, as in Nf = 2 simulations, where degener- 
ate up and down sea quarks of mass m„d are included, or fully, as in A^j = 2 + 1 calculations, where 
strange sea quarks of mass are also incorporated. 

The main motivation for studying kaon physics off and on the lattice is to test the standard 
model, to determine some of its parameters and to constrain possible new physics scenarios. From 
a lattice perspective, kaon processes fall into three broad categories. The first are processes, such 
as leptonic and semileptonic kaon decays, for which lattice QCD methods are already providing 
high precision results. The second category corresponds to processes for which lattice calculations 
are deUvering results with errors on the level of 10% or less, such as for K^-K^ mixing matrix 
elements. The last category of processes are those for which lattice calculations have failed up 
until now to provide reliable answers. Amongst them are the A/ = 1/2 rule and, more critically, 
direct CP violation in A' —> TTTT decays. 

Another motivation for studying kaons physics on the lattice is the overlap this physics has with 
chiral perturbation theory (ChPT). ChPT describes the low-energy dynamics of the pseudo-Nambu- 
Goldstone bosons of chiral symmetry breaking and has been successful in many phenomenological 
applications. Moreover, it is a very useful tool for understanding the dependence of lattice results 
on light quark masses and on volume. Recent Nf = 2 and 2 + 1 calculations, which include pions 
with masses Mji < 350 MeV, are not only using ChPT but are also beginning to provide information 
about ChPT in return. 

The talk begins with a critical discussion of the role that ChPT and other expansions can 
play in interpolating and extrapolating lattice QCD results to the physical mass point, {mi,d,ms) = 
(m^^,mf''), in view of the quark mass values currently reached in lattice calculations. In an aside, 
I present a scheme for visualizing the extent to which a lattice calculation includes the different 
effects necessary for computing a quantity of interest reliably, and a procedure for averaging lat- 
tice results. This is followed by a review of calculations of quantities relevant for leptonic and 
semileptonic kaon decays, as well as for neutral kaon mixing and CP violation in A' —> TTTT decays. 

2. Reaching the physical mass point 

Using today's algorithms, it is straightforward to perform Nf = 2 + \ calculations with a 
strange quark whose mass is around its physical value. The physical strange quark mass point 
is thus recovered simply by interpolation. 

Reaching the physical up and down quark mass point is much more difficult. Though the 
results of PACS-CS announce that calculations will soon be done directly at this point in phys- 
ically large volumes, for the moment all other simulations are being performed with larger quark 
masses. Thus, reaching the physical point still requires conducting a number of computationally 
intensive calculations with rriud < m'^^ ~ m's^ /2, extending preferably below mf'Yl2, and per- 
forming a delicate extrapolation in m„j to m^^ ~ mf'Y26. 
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To guide the interpolation to and extrapolation to m^^, a natural candidate is 5?7(3) ChPT, 
since it provides a concise framework for describing the dependence of hadronic quantities on the 
masses of the up, down and strange quarks. Moreover, ChPT in its various quenched and partially- 
quenched guises has served the lattice community well. Nevertheless, lattice calculations are reach- 
ing regions of parameter space and precisions never attained before, and it is worth considering the 
following two questions candidly: 

• What is the best way to interpolate to = m's^l 

• What is the best way to extrapolate from m^^ / 12 < rriud < rrf^ ~ m^^ /2 to rriud = 'n^l 
There are, I believe, three physically motivated options to choose from: 

(1) As already mentioned, SU (3) ChPT is a natural candidate. It has the advantage of address- 
ing both problems together, within a compact and constrained framework. Its drawback is 
that it provides similar solutions to two problems which are of a quite different nature: the 
first concerns a simple interpolation rather far away from the chiral point while the second 
involves a difficult extrapolation which reaches much deeper into the chiral regime. 

(2) SU{2) ChPT provides a means of distinguishing these two problems. For the extrapolation 
in mud, it brings to bear all of the power of chiral expansions. The interpolation in m, is not 
directly addressed, but it suffices to supplement the chiral expansion with a regular mass-or 
what I call "^avor"-expansion about m§'\ and to perform a simple polynomial interpolation. 

(3) The idea of a. flavor expansion can also be applied to the extrapolation in m„£/. To reduce 
uncertainties, this expansion should be performed about the midpoint of the interval between 
the physical point and the largest up and down quark mass considered, i.e. fhud = [m^ + 
m™j'^]/2. In this scheme, both the extrapolation in ruud and the interpolation in can be 
performed with polynomial flavor expansions. 

Let us now review these three alternatives in more detail. 

2.1 SU{3) versus SU(2) ChPT and flavor expansions: what's the difference? 

The flavor expansions are performed about regular points m„j and (i.e. they are Taylor 

expansions). This is not the case for the chiral expansions. SU (2) ChPT is an expansion about the 
singular point {mud,ms) = (0,mf''). SU (3) ChPT makes the additional assumption that the strange 
quark is chiral so that the expansion is around {mud,ms) = (0,0). 

In flavor expansions of quantities which do not vanish in the SU{2) chiral Umit, it is the 
"distance" from the expansion points, rhud or mf'', in units of the QCD scale, which determines how 
well the series converges (hence my use of the adjective "flavor"). Thus, the expansion parameters 
are (m„d — niud) /Mqcd and (mj — mf'') /Mqcd, where Mqcd ~ 1 GeV is a typical QCD scale. On 
the other hand, SU (3) ChPT expressions are expansions in mudsl with A;(; ~ 471,^^ = 0{Mqcd) 
the chiral symmetry breaking scale. In SU (2), the expansions are in mud/rus and mud/^x- 

Because ruud and are not measured directly in experiment, it is convenient to replace these 
masses by observables which are sensitive to them. ChPT suggests that M^^ and = [M^ — 
Ml/lYI^, with M^'' ~ 135MeV and M|'^'' ~ 486MeV, are particularly appropriate. Indeed, LO 
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ChPT yields = IBiriud and {M^Y = ^^s, with B = 0(Mqcd)- In terms of these variable, 
the SU (2) ChPT expansion parameters can be written {Mn/V2M^)'^ and {M^/Aj^f-, while SU{2>) 
ChPT is an expansion in {M-ji j^ -q /A^)^. Similarly, the flavor expansion parameters become A-ji = 
[Ml -Ml) /2M^cD and Ak = [{Mlf - (M|''''y]/M^cD- It is worth noting that this definition for 
remains appropriate if one assumes that itself obeys a flavor expansion in m^, i.e. = 
M^''''[l +CK{ms — ms'^)/MQCD +h.o.t], with Ck a constant. Indeed, in that case we also have 
A^ = 0[{ms — m^^)/MQCD\- On the other hand, M^j's flavor expansion in mud, = M^[l + 
2Cji{mud — rhud)/MQCD + h.o.t], is poorly behaved for the range of mud currently considered in 
lattice calculations, since the NLO plus higher order terms can be 50% or more of the LO term. 
However, this fact does not invalidate the use of flavor expansions in m„j for quantities which do 
not vanish in the SU (2) chiral limit. It merely signals that, in current calculations, the relative 
variation in is large, while the change in Mj^ with respect to Mqcd remains small. 

The expected accuracy at NLO in the SU{2) expansion around the physical mass point is much 
better than for the SU(3) case. Indeed, in SU(2) this accuracy is given by {M^^ / V2M^'''^f ~ 0. 1% 
whereas it is expected to be {M^^ /AnFjiY 5% in the SU (3) case. However, with pions of about 
450 MeV floating around, as in present day simulations, the SU (2) figure becomes {Mji/ \/2M^)'* ~ 
20%, which is much less impressive. Nevertheless, this expansion has the advantage that its con- 
vergence improves rapidly as M„ is reduced, while the SU{3) expansion parameter {My^ / AnF-^Y 
does not decrease significantly with Mj^. 

The accuracy of the. flavor interpolation in strange quark mass is generically very high. Sup- 
pose that one has performed the calculation for at least two values of the strange quark mass that 
bracket mf'' with a total spread of about 10%. The expansion parameter is then |A/f | ~ 0.01. As- 
suming that the error due to the truncation of the interpolating polynomial is on the order of the 
first omitted term, the systematic error associated with a Unear interpolation in (M^)-^ (i.e. a linear 
interpolation) will have an accuracy on the order of A| ~ 0.01%. 

In current lattice calculations, the flavor expansion in up and down quark mass is not as good. 
Assuming that we consider only pions with Mj^ < M^"^ = 450 MeV, the expansion parameter is 
|A;i-| ^ 0-05. This means that a linear extrapolation will have a truncation uncertainty on the order 
of A^~0.3% (with a coefficient that increases with u/d content). Moreover, it is straightforward to 
show that, with a quadratic flavor expansion, one can fit a chiral logarithm which gives a correction 
of up to 30% as M„ varies in the range from M^^ to M^, with a systematic accuracy better than 
0.5%. So, even in the presence of a chiral logarithm, a flavor expansion can be used. 

Let me now add a few words about the possible outcomes of implementing the different ap- 
proaches. SU (3) ChPT provides functional forms which are more constrained, i.e. which have 
less parameters, at a given order, than the SU (2) chiral and flavor expansions. That is one reason 
why SU{3) ChPT might be appealing. So let me assume, for the moment that we are fitting lattice 
results to SU(3) ChPT expressions. As M,, is lowered below \/2M| with fixed m^, SU{3) ChPT 
turns into SU (2) ChPT, except that the extended symmetry of the SU (3) theory imposes constraints 
amongst the SU (2) LECs. These constraints can be released by adding NNLO and higher terms to 
the SU (3) expansion. If theM|/A2 expansion in the SU (3) theory behaves well, then the LECs 
obtained with the fits may be SU (3) LECs of QCD, as defined in the SU (3) chiral Umit. However, 
if the assumption that the strange quark is chiral is not borne out in practice, a good fit may still be 
obtained by adding higher order terms, but the fitted LECs will most likely not be QCD's LECs. 
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In that case, one may still find that the component of the SU (3) chiral expansion is reasonably 
well behaved. If this is so, an SU (2) chiral fit ought to work and should give the SU (2) LECs of 
QCD. However, the expansion may still behave poorly for heavier pions because in that case the 
expansion parameter {M^/ \/2M^)^ may not be small. Alternatively one may use the flavor expan- 
sion approach. It deals with the strange quark mass interpolation in the same way as SU (2) ChPT, 
but differs in the choice of expansion point for the extrapolation in ChPT expands observables 
about Mji = 0, which is further from the lightest simulated than is the physical point. The 
flavor expansion, on the other hand, is performed about a value of Mj^ = which is between the 
heaviest simulated Mji and the physical value. Thus, the flavor expansion will be better behaved, 
though generically less constrained. 

ChPT is a worthy object of study in its own right, with applications which go beyond present 
lattice QCD capabilities. Thus, it is important to test its range of validity and its accuracy where it is 
applicable. It is also important to determine its LECs, since these can be used to make predictions 
in a variety of processes. However, if the goal is to determine the value of an observable at the 
physical point, one should remain agnostic in regards to the expansion used and pick the one which 
gives the lowest combined statistical and systematic error. Moreover, if the goal is to obtain the 
LECs of QCD, it may be necessary to perform calculations closer to the chiral hmit, especially in 
the case of 5f/(3) ChPT. 



2.2 SU (3) versus SU (2) ChPT and flavor expansions: examples 

To further clarify the difference between the different expansions and their applicability to 
lattice calculations, it is useful to turn to a concrete example. We consider here the expansions of 
the pion and kaon decay constants, and Fk, at NLO. In the SU (3) theory, we have [0]: 



Fn =FA\ 



FA 1 



1 



(471^3)2 
1 

(471^3)2 



+4L4(m)t| 



Xi{Ml) + \xi{Ml) 



+ 4(L5+L4)(At)-^ + 8L4(M)-^ 



+ 4(L5+2L4)(At)-f 



(2.1) 



(2.2) 



where Xn{M^) = M^"ln{M^ / jx^) and where F3 is the pion decay constant in the A^y = 3 chiral hmit. 
The up-down and strange quark mass-dependence of these two quantities are obtained here in terms 
of only three parameters: F3, L4 and L5. 

The SU (2) theory is much less frugal with parameters. At NLO it predicts [^-|5|]: 

1 



F^ = F2(l + aMjf)<! 1 
Fk = F2^{\ + a^AK 



1 



(471^2)2 
1 



[Xi{Ml)-i4{^)Ml] \ + 0{mIAk) 

"3 



(471^2)2 



;Zl{Ml)-if{^)Ml 



+ o(mX) 



(2.3) 
(2.4) 



where F2 and F2 are the pion and kaon decay constants, respectively, in the Nf = 2 chiral limit and 
where I have included a strange quark mass dependence. Thus, the SU (2) description of the mass- 
dependence of the two decay constants requires at least six parameters (F2, £4, CCp, F^ 
eight if 0{M^Ak) terms are required. 
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Figure 1: RBC/UKQCD's partially quenched pion decay constants versus pionmass squared, for two values 
of the sea pion mass, M^i = 331 ant/ 419 MeV. The only points where sea and valence quarks are degenerate 
are the square and circle with crosses. Also shown are the unitary SU{3) and SU(2) fits. Conventions are 
such that fj^ = ^/2F„ = 131 MeV. 

This number of parameters is comparable to that required in the flavor expansion of Fk and 
Fj^. Six parameters are needed if the dependence turns out to be linear and eight if curvature is 
observed, corresponding to an expansion to 0{A^,Ak). 

Let us now investigate how these considerations play out with real lattice results. I begin with 
a partially quenched, Nf = 2 + 1 study of Fji and Fk performed by RBC/UKQCD whose results 
were presented at this conference by E. Scholz [^]. These results are shown in Fig. |I|, where the 
pion decay constant is plotted against the valence pion mass squared for two values of the sea pion 
mass (331 MeV and 419 MeV). Details of the simulation are given below in Table 

In their calculation, the SU{3) ChPT expansion parameters are, at M""*^ = 419MeV: {M^"-y 
AnFtf ~ 0.1 and (M^/ AnFtfc^O.'i. The SU{2) expansion at M™"^ is not any better: (M™''^/ 
\/2M^''''')^ ~ 0.4. Thus it is not clear, a priori, which of the two expansions is better at the top of the 
range. Of course, as already mentioned, as Mj^ decreases the SU (2) expansion improves rapidly 
whereas the SU{?>) expansion parameter, (M/f ,j/47rFjf'')^ stays roughly constant. Assuming that 
SU (3) ChPT is applicable, they find very large NLO corrections to the pion decay constant, even 
at their lightest unitary point, Mj^ = 311 MeV, where they are of order 70%. They also find that the 
NLO forms do not describe their kaon results, where the down quark is replaced by a strange. This 
is perhaps not too surprising since their kaons have masses of up to approximately 570 MeV. 

With SU (2) ChPT, on the other hand, they obtain good fits and find much more reasonable 
NLO corrections, that are on the order of 30% at M^^ = 311 MeV. They use this information, 
together with that obtained from fits with partial NNLO terms and more massive pions, to conclude 
that SU (3) ChPT fails in the range of masses explored, while SU (2) ChPT is reliable. 

A few comments are in order. The first is that the fits do not take into account correlations 
which are obviously strong at fixed sea quark mass. This makes getting a meaningful figure of 
merit for the fits difficult. The second is that the results display none of the logarithmic behavior 
which becomes relevant in the extrapolation to physical Mj^: at NLO in partially quenched ChPT, 
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Figure 2: From PACS-CS, Fk as a function of the isospin averaged up and down quark masses, both 
in lattice units, at a single value of the lattice spacing, a ~ 0.09 fm /^. The black circles are the decay 
constants obtained from the different simulations, corresponding to — 156, 296, 385, 411, 570 and 
702MeV. The left panel shows these decays constants together with results from a variety ofNLO SU(3) 
ChPTfits (triangles) while the right panel displays the same results with NLO SU(2) ChPTfits (triangles). 
If all fits were good, triangles at each m^^' would sit atop the corresponding circle. In their conventions, 
fK = V2Fk = 159M6V. 



the dependence on valence quark mass is linear and with only two values of the sea quark mass, 
one cannot distinguish between a straight line and a chiral logarithm. Thus, the lattice results are 
not inconsistent with SU (2) ChPT, but they cannot be claimed, either, to confirm the relevance 
of this expansion in the quark mass range considered. Moreover, the value of Fj^ obtained by 
Unear fit would be significantly larger than the one found in the plot, though consistent within the 
final systematic error quoted by the authors. Finally, it should be remembered that the analysis is 
performed at a single, rather large value of the lattice spacing (a ~ O.llfm), and mass dependent 
discretization errors could distort the physical chiral behavior. 

PACS-CS has also investigated the applicability of the two variants of ChPT to their results 
for the decay constants [|l]], as shown at this conference by Y. Kuramashi [^. Their calculations 
are performed for six different values of the pion mass, ranging from ~ 700 MeV all the way down 
to 156 MeV. Moreover, they consider only unitary points, for which valence and sea quarks of the 
same flavor have identical masses (i.e. no partial quenching). The parameters of their simulations 
are given below in Table [l[ Their studies of the dependence of Fk on the isospin averaged up 
and down quark mass, m„rf, are shown in Fig. ^. The left hand panel displays the decay constants 
obtained directly from the simulations together with the values of these constants which result 
from fitting the simulation data to various SU (3) ChPT forms. The fits are restricted to points with 

< 410MeV. The fit results above this point are extrapolations. They find that NLO SU{3) 
ChPT fails to reproduce the dependence of Fk above Mj^ ~ 400 MeV. Moreover, they find that 
it fails to predict the strange quark mass dependence of Fk around mf'' and for Mji ~ 400 MeV. 

Again, the situation is quite different for SU (2) ChPT fits. There they find that the m„j depen- 
dence is well reproduced up to Mji ~ 4 10 MeV and only deviates from the simulation result by 5% 
at Mjc ^ 570 MeV. Moreover the wiy dependence is correctly reproduced, as it should since there 
are two mv values and this dependence is fitted by a line. 

These calculations, performed almost all the way down to the physical point, are a real prowess. 
For the moment, however, the volume considered for their lightest pion (Mji ~ 156MeV) is small. 
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corresponding to LM^ ~ 2.3. This may make it difficult to control finite- volume effects at low M^. 
Moreover, the calculations have only been performed at a single lattice spacing for now, so that 
alterations of the mass dependence by discretization errors have not yet been investigated. 

Combining the experiences of RBC/UKQCD and PACS-CS, the following conclusion seems 
to emerge: SU (3) ChPT appears to break down at the physical strange quark mass, at least in the 
presence of heavier up and down quarks, with masses larger than mud — ^^w^^ — 'wf''/3, corre- 
sponding to > 400 MeV. 

3. Aside on a classification sclieme for lattice simulations and on the averaging of 
lattice results 

Before turning to quantities of phenomenological interest, I wish to say a few words about the 
methodology that I will follow in reviewing lattice results. 

3.1 Of stars and lattice calculations 

The FLAVlAnet Lattice Averaging Group (FLAG) is in the process of putting together a clas- 
sification scheme for lattice calculations. The goal is to provide tables which, at a glance, give 
the reader a sense of how thoroughly a given calculation includes all of the necessary ingredients, 
based on a list of pre-defined, objective criteria. Since this collective work has not yet been final- 
ized, 1 propose a personal version of the scheme here. It is based on a starring system, reminiscent 
of the one used in a famous, red restaurant guide: 

-k-k-k indicates that this aspect of the calculation is fully satisfactory; 

indicates that the corresponding ingredient has not been fully included, but that the invstiga- 
tions performed allow for a reasonable estimate of the ensuing systematic error; 

• indicates that the calculations performed are not sufficient to provide a reliable estimate of 
what is missing. 

More specifically, here are the criteria which I use for starring the calculations reviewed below: 

• publication status 

published 
irk preprint 
★ proceedings, talk 

• action, unitarity 

•k-k-k local action, unitary calculation 
irk non-local action and/or discretization induced unitarity violations 

• flavors 

all dynamical flavors required for the process under study are included 
irk some dynamical flavors missing, but at least Nf>2 
i< Nf = (i.e. quenched calculation) 

• renormalization 

nonperturbative with nonperturbative running 
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nonperturbative with perturbative running at GeV energies, or perturbative at two-loops 
or more 

•k one-loop perturbative and/or discretization which leads to poorly controlled operator 

mixing 

• extrapolation/interpolation to physical mass point 

Let M'^'iyp be a mass that is representative ("typical") of the masses of the hghtest pion 
variants that contribute to the mud dependence of the quantities studied. ^ 

M^^Pyp < 200MeV with NLO or better ChPT or any other demonstrably controlled 

functional mass dependence 

M^'Pyp < 350 MeV and reliable estimate of extrapolation error 

* M™^p > 350MeV 

• continuum extrapolation 

•k-k-k 3 or more lattice spacings with at least one a < 0.08 fm and controlled scaling 
2 or more lattice spacings with one a < 0. 1 fm 

* a single lattice spacing or all a > 0. 1 fm 

• finite volume 

■k-k-k LM^'" > 4 (and numerical volume scaling study) 
-kk 3 < LM^'" < 4 and well motived analytical corrections 

* LM^"^ < 3 or 3 < LM^"^ < 4 and no quantification of finite-volume effects 

where M™" is the mass of the lightest pion contributing to finite-volume effects. 

3.2 Averaging of lattice results 

Now that results for various quantities of phenomenological interest are emerging from lattice 
calculations in which most effects are realistically taken into account, it is important to set forth 
objective, quantity independent averaging procedures. In particular, that means taking literally the 
statistical and systematic error estimates provided by the authors in a refereed pubUcation. It also 
means only considering calculations in which all relevant sources of systematic uncertainty have 
been accounted for. Since we are still in the early days of reaUstic lattice calculations, this rule 
might have to be bent slightly at first to include results which are close to reaching this goal. 

The averaging procedure which I adopt is the following. Given a list of results which sat- 
isfy the basic criteria described above, I perform their weighted average, with an inverse weight 
obtained by adding the statistical and systematic covariance matrices in quadrature. To determine 
the statistcal error on the average, I construct a X'^ with only the statistical correlation matrix and 
perform a standard Ax^ analysis. For the systematic error, since one does not generically expect 
them to compensate from one calculation to the next, I take the smallest total systematic uncer- 
tainty amongst those obtained in the most complete calculations. In cases where either statistical 
or systematic errors are not symmetric, I symmetrize them. 

' This "typical" mass depends on the fermion formulation used, on the quantities studied, etc. Since it is meant 
to be indicative, I have kept its determination simple. For staggered calculations I take the RMS of the masses of 
the different tastes; for non-staggered on staggered, the RMS of the valence and the sea taste-singlet pion masses; for 
Wilson, overlap, domain-wall, the RMS of the valence and sea pion masses (i.e. simply the lightest pion mass for 
unitary calculations); and for tmQCD, I have taken the charged pion mass, though some sort of isospin averaging should 
probably be performed. I thank C. Aubin, J. Laiho, S. Sharpe and R. Van de Water for enlightening correspondence. 
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There will be some statistical correlations between results obtained from the same set or from 
subsets of a given ensemble of gauge configurations. There will also be some correlations in the 
systematic errors of calculations which make use of similar methods. However, such correlations 
have not yet been analyzed in any detail and I choose to neglect them here when computing world 
averages. For computing an average's statistical error, though, I keep only the statistical error of 
the calculation, amongst those performed on a same set of configurations, that makes use of the 
largest fraction of these configurations. Correlations may be added more systematically later, once 
they are better understood. 

In situations where some results have significantly smaller systematic uncertainties, for reasons 
which are not fully understood, one can provide an average with and without those results. 



4. \Vus\ from experiment and the lattice 

A precise determination of the magnitude of the CKM matrix element Vus allows for a precision 
test of CKM unitarity as well as of quark-lepton universality and provides constraints on new 
physics, through: 



G 



1 + 0(^1, (4.1) 



where iyud^Vus^Vub) forms the first row of the CKM matrix and where Gq is the Fermi constant 
as measured in quark decays, whereas G^ is the same constant as determined in muon decays. 
Eq. ( |4.l| ) accounts for the fact that what is actually measured are not the CKM matrix elements, 
Vqqi , but (Gy/G^ ) X I Vqqi \ . The large amounts of new experimental results from BNL-E685, KLOE, 
KTEV, ISTRA+ and NA48 provide the opportunity for testing this aspect of the standard model 
with unprecedented accuracy. 

The current situation on the measurement of the relevant CKM matrix elements is: 

• \Vud\ = 0.97425(22) [0.02%] from nuclear jS decays ||] 

• \Vus\ = 0.2246(12) [0.5%] from Kb decays |g] 

• \Vus/Vud\ =0.2321(15) [0.6%] from i^f2 decays [g 



Vub\ = 3.87(47) • 10"-^ [12%] from exclusive and inclusive b uiv decays ^ 



where a factor of {Gq/ G^) is implicit, as per Eq. ( |4. 1[ ), and where the percentages in square brackets 



indicate, for convenience, the relative error of the measurement. 

The Flavianet Kaon Working Group combined the first three measurement to squeeze out a 
Uttle additional precision on [Vusl I have updated their analysis here to take into account the 
new result for \Vub\ 

• \Vud \ = 0.97425(22) [0.02%], which implies the following contribution to the uncertainty in 
Eq.(0): 5|Kd|' = 4.3 -10-4, 

• \Vus\ = 0.2252(9) [0.4%], which impUes the following contribution to the uncertainty in 
Eq.(y): 5|K,vP= 4.2 -10-4, 

• and the contribution from Vi,b to Eq. (1.1), |V„tp ~ 1.5 • 10^^, is so small that its error bar is 



irrelevant. 
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Ref. Nf action 


a[fm] 


LM'li'" M™'"[MeV] 
typ/val typ/val 




PDG'08 [11] 






1.193(6) 



ETM'08 [12] 



r||LQCD'06 
MILC104-07 

[Bill 

HPQCD/'Q7 
UKQCD ^ 

RBC/'08 ^ 
UKQCD [H 

ALV'08 

PACS-CS'08|[Ill 

BMW 08 



2 tmQCD 0.07, 0.09, 0.10[F;;r] 
0.13[ro] 

AsqTad 0.06,0.09,0.12, 



2+1 DWF/ 

J>-^MILC 



2+1 KS^iLC o_i5[^^] 
2+lKS™LC 0.09,0.12,0.15[T] 
0.11[n] 

0.09,0.12[T/F;,] 

0.09 [a] 
0.065,0.085, 
0.125[S] 



2+1 DWF 
2+1 DWF/ 
2+1 NP-SW 
2+1 SW 



3.6/3.6 


260/260 


1.196(1 3)(7)(8) 


5.1/3.5 


420/290 


1.218(2)+^^ 


5.3/4.2 


300/240 


1.197(3)+f3 


4.8/4.1 


360/310 


1.189(2)(7) 


4.1/3.4 


290/240 


1.205(18X62) 


5.3/4.2 


300/240 


1.191(16)(17) 


2.3/2.3 


160/160 


1.189(20) 


4/4 


190/190 


1.19(1X1) 



Table 1: Parameters of the simulations used by various collaborations for calculating Fk /F-,1, together with 
their result for that quantity ( results in italic were presented at this conference ). The column Nf indicates the 
number of sea quark flavors considered. The symbols in brackets in the a[fm] column indicate the quantity 
used to set the scale. Also given are the masses of the lightest valence pion simulated, as well as the "typical" 
lightest pion mass defined in footnote^ 

At the time of the conference, |V„.s| was no longer the dominant source of uncertainty in Eq. (p~l|). 
However, the new result for \Vud\ [pll makes it a dead heat. Combining these results yields: 



\Vud\^ + \Vus? + \Vub\^ = 0.9999(6) [0.06%] . (4.2) 



This result is fully consistent with the standard model. However, within one standard deviation, 
new physics at a scale Kj^p > 3 TeV cannot be excluded and within three standard deviations, this 
scale drops down to Aatp > 2 TeV. 

4.1 \Vus/Vud\ from K,7i^ jxv 

In 2004, Marciano pointed out a window of opportunity for determining iKs/Kt/l from the 
ratio of leptonic decay rates r{K /xv(7))/r(7r Mv(7)) [|l^. Calculating 0(a) radiative cor- 



rections to this ratio, he obtained (see update in fll IQ): 

It/ I ^ 

j-^-^ =0.2757(7) [0.25%] . (4.3) 

\Vi,d\ Fji 

Thus, a precise lattice calculation of Fk/Fji will allow a high precision determination of |Vi„ /Vurf|- 
One needs to determine Fk/Fj^ to: 

• 0.5% to match the uncertainty on \Vus\ obtained mK ^ niv decays, 

• 0.25% to match the experimental uncertainty in r{K /xv(7))/r(7r Mv(7)). 
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Ref. 

ETM'08 



NPLQCD'06 [13] 



MILC 04-07 1 14, 15] 



HP/UKQCD'07 [[T§] 
RBC/UKQCD'08 ^ 
ALV'08 [jl8p 
PACS-CS'OB ^ 
BMW'08 



irk 

-k 



irk 
irk 
irirk: 
irk 

kkir 



kkr 
kkk 
kkk 
kkk 
kkk 
kirk 
kkk 
kkk 



irk 
k 
irk 
k 
irk 
irk 
kkk 
kkk 



irk 
k 
kkk 



irk 
kkk 
kkk 
kkk 
kkk 
kkk 

k 
kkk 



Table 2: Starring of the simulations used to obtain Fk/Fji, according to the criteria put forth in Sec. 3.1 



Pk/Ftz is an SU (3)-flavor breaking effect, i.e. 

(mI-mI\ 

Fk/F^ = \ + 0[ M (4.4) 

and it is the deviation from unity that we are actually calculating, which makes the target accuracies 
a little less forbidding. 

In Table |I| I summarize the parameters and results of all unquenched lattice calculation of 
Fk/Fti- The corresponding consumer report is given in Table 



Of all these calculations, the most advanced is that of MILC []14|, |150, but the calculation of the 
BMW collaboration, presented at this conference by S. Diirr, should rival it once completed. The 
calculation of PACS-CS performed very close to the physical up and down quark mass holds 
great promise. However, as it stands, it is missing a continuum extrapolation and may also suffer 
from significant finite-volume errors. 

To illustrate lattice results for Fk/Fj^, let me briefly present those of BMW. The calculations 
are performed in volumes up to 4 fm, with pions as light as 190 MeV and lattice spacings down to 
0.065 fm. The parameters of the calculation are summarized in Table |l|, and details of the ensembles 
can be found in [p^. The results are plotted in Fig. ^, as a function of M| in physical units, with 
the scale set by the E mass as in [[l9|]. The plot shows the extrapolation of the results for Fk/Fji 
in to the physical point. A large variety of functional forms have been tried, ranging from 
NLO SU (2) ChPT to polynomial expansions. Three different cuts on pion mass have been made: 

< 420 MeV, 470 MeV and 600 MeV. The continuum and mass extrapolations are combined. 



by allowing for the parameters of the functional mass dependence to acquire or a corrections. 



Finite- volume effects are subtracted at two-loops in ChPT, using the results of [20]. The procedure 
for estimating statistical and systematic uncertainties is very similar to that in [19]. It should be 
noted that the shift in Fk/Fji from the lightest pion mass to the physical point is less than 2%. The 
preliminary result is given in Table [T[ 

Unquenched, lattice results for Fk/F^ are summarized in Fig. ^, where my average for this 
quantity, obtained as explained in Sec. is also given. This average includes only the published 
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1.25^ 
1.2- 

1.1- 
1.05- 



a«0.125 fm 
a=0.085 fm 
a=0.065 fm 
a=0 




0.2 

M^[GeV^] 



FigureS: Fk/F-^ versusM^ in physical units, as presented by BMW at this conference. The different symbols 
correspond to different lattice spacings, as indicated in the plot legend. The curves represent the result of a 
combined chiral and continuum extrapolation fit for each lattice spacing, as well as in the continuum limit. 
The results have already been interpolated in strange quark mass to the physical value. The particular fit 
shown corresponds to a NLO, SU (3) ChPTfit with < 470 MeV and a^ discretization errors. 



Nj = 2+1 
(MILC) 



T ' r 

ETMC '07 



NPLQCD '06 (•) 

HPQCD/ 
UKQCD '07 (•) 

MILC '07 (») 

ALV '08 

RBC/UKQCD '08 (•) 
PACS-CS '08 
BMW 'OS 



F^/F^=1.194(3)(10) 

J , I , 



Figure 4: Summary of unquenched lattice results for Fk/Fji, together with my average. The latter is 



obtained as described in Sec. 3.2 and in the text. The results marked with a "(*)" are those included in the 



average. The smallest error bar on each point is the statistical error and the larger one, the statistical and 
systematic errors combined in quadrature. The references are as in Table^ 



Nf = 2+\ results [13-16, ^] in which many systematic uncertainties have been estimated. The 
systematic error is taken from [15]. The total uncertainty on this quantity is 8{Fk/Fji)^^" = 0.8%. 
This corresponds to an uncertainty of d{FK/Fji — l)'"' ~ 5% on the calculated SU (3)-flavor break- 
ing effect, which is much better than the accuracy obtained on the SU (3)-flavor breaking in the 
form factor for K —>■ Tiiv, 5/+(0) ~ 15%. Nevertheless, this uncertainty still leads to a larger the- 
ory error in the determination of \Vus\, i.e. 0.8% vs 0.5%. Since Fk/Fji is a straightfoward quantity 
to calculate, one may expect steady improvements in its lattice determination, especially in light of 
the recent progress by PACS-CS 
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4.2 \Vus\ ^rom K^Tiiv 

K niv decays provide an alternative way to determine \Vi,s\. This measurement requires the 
theoretical calculation of the vector form factor, f+{q^), defined through: 

{n+{p')\ur^s\K''{p)) = {p + p'-q ^ , " U+(^')+^M " 2 " /o(^') ' (4-5) 

with q = p — p'. The best precision is currently obtained by measuring the form factor shape in 
experiment and extracting, from the total rate [^, 

\Vus\ X /+(0) = 0.21664(48) [0.22%] . (4.6) 

The experimental error is of similar size as in the ratio of leptonic kaon to pion decay rates. To 
fully exploit the experimental results requires a determination of /+(0) to 0.22%. 

To the extent that it is applicable here, the theoretical framework for attacking this problem is 
SU{3) ChPT [|, |l|]. The chiral expansion for /+(0) is given by: 

/+(0) = l+/2 + /4 + --- , (4.7) 

where the AdemoUo-Gatto theorem [ ^2| ] and ChPT yield: 

Thus, there are no contributions for the 0{p'^) LECs and this NLO contribution is fully determined 
by Mk, and Fj,. 

This means that a sub-percent level determination of /+ (0) requires a calculation of NNLO 
and higher corrections, since 

A/^/+(0) - 1 -/2 = o( ^^^\f'^^ ^ -3% (4.9) 

is comparable in size to /2. To fully exploit the experimental accuary "only" requires an accuracy 
of 7% in the calculation of A/. 

What is known about and more generally A/? The NNLO chiral logs have been com- 
puted [^, ^], and they require 0{p^) LECs for renormalization. Estimates have been made 
of these LECs [^-27] and in p4[ ] it is shown that they can be determined from the slope and 
curvature of f+{q^). The reference value for A/ is still taken to be the quark model result, 
A/ = -0.016(8) lllll- 

In Table |, I summarize the parameters and results of all unquenched lattice calculations of 
/+(0). The corresponding consumer report is given in Table Q 

The lattice methodology for the calculation of /+(0) — 1 was set forth in [^]. It consists of 
three main steps: 

1. Use a double ratio of three-point functions to obtain: 

2 2VMm~n{n\Vo\K){K%\K) 
mq,na.) Mk + M, {K\Vo\n){K\Vo\K) ' ^ " ^ 

This yields a determination of foiqliax) with a statistical error less than about 0.1%! 
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Ref. 



N 



f 



action <3[fm] L[fm] 



typ/val 



/+(0) 



JLQCD'05 [28] 



RBC'06 [29] 



ETM'08 [|30|] 
FNAL/MILC'04 [||] 



RBC/UKQCD'07|32] 



2 
2 
2 

2+1 
2+1 



NP SW 

DWF 

tmQCD 

KS+Wil 

DWF 



0.09 
0.12 
0.11 



1.8 
2.5 
2.7 



0.11 1.8,2.8 



550/550 0.967(6) 
490/490 0.968(9)(6) 
260/260 0.957(5) 

0.962(6)(9) 
290/240 0.9644(33)(34)(14) 



Table 3: Parameters of the simulations used by various collaborations for calculating f+(0), together with 
their result for that quantity. The description of the columns can be inferred from the one given in Table ^ 



Ref. 








JLQCD'05 Q 
RBC'06 [29] 
ETM'08 |3C] 
FNAL/MILC'04 [31] 
RBC/UKQCD'07[||] 









i^-k-k 



Table 4: Starring of the simulations used to obtain /+(0), according to the criteria put forth in Sec. 3J_ 



2. Compute fo{q^) at various and use an ansatz to interpolate and get /+(0) = /o(0). 

3. Interpolate/extrapolate in light quark mass to the physical mass point. 



RBC/UKQCD | |32| ] actually combine steps 2 and 3, using the functional form: 

1 +f2{MKMn) + {Ml-Mlf{Ao+A,{Ml + Ml)) 



fo{q\MK,M^) 



l-q2/{Ma+M,{Ml+Ml)Y 



(4.11) 



where Aq, A i , Mq and M\ are parameters and where a polynomial ansatz was made for NNLO terms. 
This combined fit is shown in the two panels of Fig. |5[ Their results fit 1 + f2{MK,Mj^) + NNLO 
well, though their fits do not take correlations into account. The claim that they are sensitive 
to NNLO effects seems to be justified. Moreover, the extrapolated result is only two standard 
deviations below the result obtained at their lightest pion mass and the claimed error on /+ (0) — 1 
is a rather conservative 14%. The caveats are that is approximately 15% too high and the 
calculations were performed at a single, rather coarse lattice spacing of a = 0.1 14(2) fm, meaning 
that discretizations errors can only be guessed. Nevertheless, this is the first convincing lattice 
calculation of /+ (0) — 1 . 

Lattice and non-lattice results for /+(0) are summarized in Fig. ^ together with the "average" 
which I obtain by copying the result of p2|]. The total uncertainty on /+(0) is 5/^(0)'"' = 0.5%. 
This means that K —>■ niv decays still give the best accuracy for \Vus\. Moreover, one can anticipate 
that the current error, 5(/+(0) -l)'"' = 14% will be reduced thanks to the use of: stochastic 
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Figure 5: Combined and chiralfit of foiq^) to Eq. fgjTp by RBC/UKQCD '08 /|2|/. The left panel 
displays the lattice values of fo{q^ ) vs q^, together with the fit curve obtained at the physical values ofMjc 
andMK- The lattice points were shifted in pion and kaon mass at fixed q^ using the fit result. The right panel 
displays the extrapolation o//+(0) = /o(0) in to physical pion mass. 



Leulwyler & Roos "84 
Bijiiens & Takivera '03 
Jarain el al '05 
Cirigliano et al '05 



SPQcdR '04 
JLQCD '05 
RBC '06 



FNAL/MILC '04 
RBC/UKQCD '07 (•) 

ff' (0) = 0.964(3)(4) 



Figure 6: Summary of lattice results for /+(0), together with the results obtained in various models. Also 



given is my average of the unquenched lattice resutls. The latter is obtained as described in Sec. 3_2 and in 
the text. The results marked with a "(*)" are those included in the average. The smallest error bar on each 
point is the statistical error and the larger one, the statistical and systematic errors combined in quadrature. 
The references for the lattice works are as in Table ^ with, in addition, SPQcdR '04 /{33|/. The others are 
Leutwyler and Roos '84 /pi]/, Bijnens and Talavera '03 /p4|/, Jamin et al '05 /|25|/, Cirigliano et al '05 /||/. 



sources, as used in partially twisted boundary conditions 

in [34, 38-40], which enable to determine f+{q^) directly at 



1^], applied to form factors 
0[p 



351. 



5. K ^ nn decays on the lattice 

The phenomenology of ^ ^ nn decays is extremely rich, and has been highly instrumental 
in developing the standard model. In the isospin limit, the amplitudes for these decays can be 
decomposed in terms of amplitudes Aje'^' , I = 0,2, where / is the isospin of the final two-pion 
state and 5/ is the strong scattering phase in that channel. CP violation implies that A| / Aj. CP 
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violation occurs in two ways in Ki decays. Ki is mostly CP odd, and decays predominently into 
three pions. But it has a small CP even component, through which it can also decay into two pions. 
This process is known as indirect CP violation, and is parametrized by 

^ r[^L^(7r7r)/=o] ,-0, . . /ImMi2^,^ 

— ^ ' J ~ e^<' sm (pe I h C, 



AMk 



(5.1) 



T[Ks ^ iKK),=o] 

with ^ = ImAo/ReAo, AM/f = M^^^ —^Ks and M12 to be defined below. decays can violate 
CP through another channel, by having its CP odd component decay directly into two pions. This 



process is known as direct CP violation, and is parametrized by [41]: 



1 T[K,^{%%)i=2] 



V2T[Ks 
1 

V2 



(7r7r)/=o] 



T[Kl ^ {kk)i=2] T[Kl ^ (7r7r)/=o] 



1 ^i(K/2+S2-So) Re^2 

ReAn 



T[Ks- 
ImA2 



{7171)1=2 
ImAo 



T[Ks ^ {7l7l)i=o] 



ReA2 ReA() 



(5.2) 



Experimentally a lot is known about these different processes [|1 IQ. The Kl-Ks mass difference 
is measured to high precision, i.e. AMk = (3.483 ±0.006) x IQ-^^MeV [0.2%]. K ^ nn decays 
exhibit a strong enhancement of the 7 = channel over the 1 = 2 channel, \Aq /A2 \ ^ 22.2, known as 
the A/ = 1/2 rule, which is still in need of an explanation after over forty years. The parameter for 
indirect CP violation has also been measured to high accuracy, |£| = (2.229±0.012)-10-3 [0.5%], 
with a phase 0e = (43.5 ± 0.7)'^ [1.6%]. And after an experimental effort of nearly thirty years, 
direct CP violation was also measured, yielding Re(£7£) = (1.65 ±0.26) • 10"^ [16%]. 

5.1 K^-EP mixing in the standard model and Bk 

K^-K^ mixing is responsible for the Kl-Ks mass difference as well as for indirect CP violation 
in K ^ Tin. In the standard model, the CP violating contribution occurs through a local AS = 2, 
four-quark operator, once the heavy, standard model degrees of freedom are integrated out. The 
corresponding amplitude is 

2MkM\2 = (^"[jr.f ='[/:") =c™(At)(^|Oi(Ai)|^*^) , 

where is a short-distance Wilson coefficient and where 

{sd)v-A{sd)v-A and (^VilAt)!^") = y^MbM ■ 



(5.3) 



In terms of theses quantities, a revised [^, 43] standard model analysis [44] yields: 
|£| ~ K,CeBK [lm(A;2)T]„5o(x,) ±2Im(A;A;)T],,5o(x,,x,) ± lm{Xf)T^ccSoixc)] , 



(5.4) 



where Cg is determined by well measured quantities, Bk = Cf^{lJ.)BK{lJ.) is the renormalization- 
group invariant B-parameter, Xq = VqdV*^ and rj^qi, So are short-distance quantities. K"g parametrizes 
the corrections to the standard analysis [0], which arise from \/2sin0e and ^ in Eq. (5J^) at leading 
order. This approximation is not necessary and should not be used for precision tests of the SM 
once ^ is known better. A rough estimate yields = 0.92(2) [|43[], which implies an 8 ± 2% 
downward shift in the SM prediction for £. From Eq. ( |5.4[ ), it is clear that a measurement of |£| and 
a determination of Bk imposes contraints on ImA*^, ImA*^ and ImA*A*, as shown in Fig. 0. 
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Given how accurately |£| is measured, one may wonder why the constraint that it gives on 
the summit of the triangle is not any better. To help answer this question, in Fig. ^ I display SM 
predictions for |e| obtained in different ways. The starting point is a global CKM fit using Eq. (5.4), 
in which the experimental measurement for |£| is not included and where B/f = 0.723(11)(35) [5%] 
(from Fig. g) and \Vcb\ = 0.04059(38) (58) [1.7%] The topmost theoretical prediction for |£| 
is obtained from this global fit, allowing all quantities to fluctuate within their error bars. ^ The 
next result is obtained by freezing \Vcb\ to its central value. The third value results from fixing Bk 
to its central value. The fourth is obtained by freezing both \Vcb\ and Bk, and the fifth by fixing 
Bk to its central value and the four CKM parameters to their best global fit values. The last is the 
experimental measurement quoted above. 

As the second point indicates, a determination of \Vch\ to infinite accuracy only reduces the 
uncertainty on the prediction for |e| from 10% to 9%. Significantly improving the accuracy on Bk 
has a similar effect, since the uncertainty on |e| is also 9% in that case. The fourth point indicates 
that the uncertainty coming from sources other than Bk and |Vcf,| is a little less than 8%. It is only 
when the uncertainties on Bk and CKM parameters are assumed to be zero that the error on the SM 
prediction for |£| falls to 5% . The latter is due to perturbative uncertainties and to the error on JCg. 

It is interesting to note that the SM prediction for |e| is now about 18% below the experimental 
value, down from what it was until recently. The decrease is mainly due to the fact that the central 
value for Bk has dropped by more than 15% over the last decade (see discussion below and caveats) 
and to the presence of the correction related to in Eq. ( 5.4). This potential tension between theory 
and observation has led the authors of [15, 43, 46] to investigate new CP violating contributions 
to = 2 observables. For the moment, as the top point in Fig. ^ indicates, the discrepancy is 
about two standard deviations. However, if uncertainties on Bk and on CKM parameters improve 
substantially, this discrepancy could become significant, as the fifth point of Fig. || indicates. 

Bk is a quantity which has a long history on and off the lattice. However, because of space con- 
straints, in Table ^ I only summarize the parameters and results of unquenched lattice calculations. 
The corresponding consumer report is given in Table |[ 



^Systematic errors due to theory are assumed to have a gaussian distribution, so that the deviation of the predicted 
from the measured lei can be counted in standard deviations. 
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1 


1 o-i+O^O 


1 

(global fit) 






(IV^J fixed) 




^■<: 


(B^ fixed) 


1 — • — 1 




(B^&IVJfixed) 




^■<: 


(B^ & CKM fixed) 


• 

1 


2.229+0.012 (expt.) 

.1.1 



1.5 



2.5 



3.5 



lei X 10 



Figure 8: Standard model predictions for |e|, using results from global CKM fits by CKMfitter / [J^ 
compared to experiment. The various predictions are obtained by fixing different inputs to their central 
values, as described in the plot legend and in the text. 



Ref. 



action 



a[fm] L[fm] 



M™'"[MeV] 
typ/val 



Bk 



JLQCD'08 

ETM'08 [|9(] 

HPQCD/ „ 
UKQCD'06 [|^] 

RBC/UKD 
'07-08 pi 



Baeetar08 |52] 



2 
2 

2+1 

2+1 
2+1 



Overlap 0.12 1.9 

OS/tmQCD 0.07,0.09 2.1,2.7 

0.125 2.5 



If c HYP 
^■^MILC 



DWF 

ifoHYP 



0.11 

> 0.06 



1.8,2.8 
4 



290/290 
300/300 

460/360 

290/240 
300/240 



0.734(5)(55) 
0.78(3) 

0.85(2)(18) 

0.717(14X35) 
5Bk ^3% 



Table 5: Parameters of the simulations used by various collaborations for calculating the renormalization 
group invariant parameter Bk, together with their result for that quantity. The values ofBK in the table have 
been obtained at NLO from the RI/MOM or MS-NDR values given in the papers. The description of the 
columns can be inferred from the one given in Table 0. 

At this conference, A. Vladikas presented new, preliminary results for Bk from ETM [p^]. 
Their calculation makes use of Osterwalder-Seiler valence quarks on the ETM, Nf = 2, tmQCD 
seas to ensure automatic 0(a) -improvement as well as multiplicative renormalization of the AS = 2, 
four-quark operator. This implies that their calculation suffers from 0{a^) unitarity violations, 
which must be controlled. Their plan is to use ETM's three lattice spacings, a ~ 0.07, 0.09, 
0.010 fm to extrapolate to the continuum limit. For the moment, though, all results are obtained 
from simulations performed at a ~ 0.09 fm. They extrapolate in up and down quark mass using 
NLO, partially-quenched SU (2) ChPT ^] and interpolate in valence strange quark mass linearly. 
Note that the extrapolated value of Bk is only ~ 3% below its value at the lightest up and down 
quark mass, suggesting that this extrapolation is well controlled. The renormalization is performed 



nonperturbatively in the RI/MOM scheme []530. The continuum limit and finite- volume corrections 
are still missing. Discretization induced unitarity violations will have to be investigated. 

Unquenched, lattice results for Bk are summarized in Fig. ^ where my average for this quan- 
tity, obtained as explained in Sec. p.2[ , is also given. Only the two Nf = 2 + 1 results [50, ^ ^ 
are taken into account in this average. The systematic error is taken from plj , ||]. However, each 
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Ref. 



0- 



JLQCD'08 [48] 



ETM'OB [|49|] 

HPQCD/ _ 
UKQCD'06 ||§ 

RBC/UKQCP 
'07-08 lUflll 



-k-k-k 



-k-k-k 



irk 



-k* 
-ki^-k 

-kir-k 



-k* 



Table 6: Starring of the simulations used to obtain Bk, according to the criteria put forth in Sec. 3.1 



JLQCD I 



HPQCD/UKQCD '06 {•) 
RBC/UKQCD'08(") 



= 0.723(11)(35) 



Figure 9: Summary ofunquenched lattice results for the renormalization group invariant Bk, together with 



my average. The latter is obtained as described in Sec. 3.2 and in the text. The smallest error bar on each 
point is the statistical error and the larger one, the statistical and systematic errors combined in quadrature. 
The results marked with a "{*) " are those included in the average. The references are as in Table | 



calculation was performed at a single, rather coarse value of the lattice spacing. This means that 
these results, and thus the average, suffer from a poorly controlled discretization errors. 

As noted above, the value of Bk has come down quite significantly compared to JLQCD's 
standard quenched value of a decade ago [Q]. In particular, (fi/f)«Bc^^V(^/f)/Lecz) ~ 0.83(8). 
This drop cannot really be ascribed to the inclusion of sea quark effects, since comparably low 
results were obtained in the quenched approximation at comparable lattice spacings [55]. However 
there, a continuum extrapolation based on two calculations performed at a ~ O.lOfm and 0.067 fm 
increased the result to (Ba:)/^^^ /^^J^)jlqcd ~ 0-90(9). Thus, it is very important to clarify this 
situation by investigating the continuum limit of B/f in 2 + 1 flavor simulations. 

The total lattice error on Bk is 5B^^* =5%, which is comparable to the other uncertainties in 
the standard model prediction for e. As discussed above, to improve this prediction and possibly 
reveal new physics, we must not only reduce the error on Bk, but also improve the determination 
of CKM parameters, of the correction related to JCg and eventually of the short distance QCD 
coefficients. 
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6. Conclusion 

Lattice QCD simulations have made tremendous progress in the last few years. 2 + 1 flavor 
lattice calculations with pion masses as low as Mji ~ 190MeV in (4fm)^ volumes, and lattice spac- 



ings down to ~ 0.065 fm have already been performed [|19p. Moreover, as PACS-CS has shown ||1[], 
simulations at physical are around the corner. Thus, it is has now become possible to reach the 
physical QCD point {Mji ~ 135 MeV, a ^ 0, L ^ oo) in a controlled fashion. 

Quantities such as Fk/F-,^ and ff°^ (0) are already being computed with percent or better 
accuracy and are having an important impact on SM and BSM tests. Quantities such as B/f are 
reaching the sub 10% accuracy level and have errors which match those from other sources. Calcu- 
lations of e' /e and the A/ = 1/2 enhancement still have 100% uncertainties despite the impressive 



Nf = 2+1 RBC/UKQCD effort [g^, but perhaps not for long Q. Many quantities are still 
missing continuum extrapolations. 

NLO SU (3) ChPT appear to be having trouble at the physical strange quark mass, at least in 
the presence of heavier up and down quarks, whereas SU (2) ChPT performs better. However, these 
inferences require further investigation, in particular once continuum limits have been investigated. 

Concerning the extrapolations and interpolations required to reach the physical mass point 
(m„j,mi) = my advice is to keep an open mind and to pick the approach which gives 

the lowest combined statistical and systematic error. 

To conclude, the age of precision, nonperturbative QCD calculations is dawning, and the next 
few years should bring many exciting results. 
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